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Critical exponents in the CP^ -1 model, which describes localized- moment ferro- and antifer- 
romagnets (N — 2 in the Heisenberg model), are calculated from two-particle Green's functions to 
first order in 1/N. For al = 2 + e the results agree with earlier renormalization group calculations. 
For d = 3 the leading 1/iV-corrections turn out to be very large at N = 2. For d — 4 — e the 1/N 
-corrections are small at any N and insufficient to describe correctly the magnetic phase transition. 
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0^ ' Recently the large- N approach has been extensively applied to investigate thermodynamic properties of localized- 
moment magnets [[y. Two ways of constructing the large- N approach are possible: in the 0(N + 1) (which is a 

,— ' ■ generalization of the nonlinear sigma model to the case N > 2) and in the SU (iV)-model. At N = 2 these models 
describe the usual Heisenberg magnet and coincide because of existence of homomorphism of SU{2) on 0(3) group: 

n = z'az (1) 

where z is the two-component complex field with the constraint condition 

|z| 2 = l (2) 



(N 



1 n is the corresponding element of 0(3) group, a are the Pauli matrices. At the same time, for N > 2 such an 
C^) equivalence does not hold, and two models under consideration correspond to two different pictures of the localized- 
■ magnet excitation spectrum. The N — > oo limit in the O(N) model corresponds to the quantum spherical model 



X 



which appeared to be rather successful for a 
is based on the SU(N) rather than 0(N) 



||, and the 1/N corrections are small Q. However the approach 
wide range of problems concerning low-dimensional magnets (see e.g 

model. The N — > co limit in the SU (iV)-model corresponds to the mean-field approximation (the modified spin- wave 
theory) [gJ.Tbis approximation was successfully applied to two-dimensional Heisenberg model |5|,|6| and lead to the 
result for correlation length, which was consistent with the one-loop scaling analysis jlO|]. Corrections of the first 
order in 1/N modify only preexponential factor in this case. But when applied to the Heisenberg model with 
d > 2, where a finite phase transition temperature exists, the mean-field approximation leads to some difficulties [§. 
In particular, it yields values of critical exponents /3 = 1, 77 = 4 — d, which are far (for not small d — 2) from the values 
obtained from the 1/N expansion in the 0(N) model M. 
Q ' In the present paper we investigate 1/N corrections to the above-mentioned critical exponents of the mean- field 
. approximation in the SU (JV) model describing quantum magnets. In the continuum limit this model reduces to 
the CP 1 ^^ 1 model The latter model enables one to use standard field theory methods for treating static and 

dynamical properties in the vicinity of the magnetic phase transition. 

To calculate the whole set of critical exponents we have to consider the non-uniform longitudinal susceptibility 
which may be represented in terms of a two-particle Green's function. The one- and two-particle Green's functions 
can be obtained from the generating functional for the CP N ~ X model JT3| 



Z[h,ht]= J z m ]5(\z\ 2 - l)exp(-5[^,2; m ]) 

l/T 

S[zl, z m ] = 2 Ps J dr J d d r [|cV| 2 - \z^d a z\ 2 + tfz + hz*] (3) 


where p s is the spin stiffness, the magnon velocity is put to be unity, z m and h m are iV-component complex fields: 
to = 1...N . Representing the delta function in (^) as a path integral over the slave field A and decoupling quartic terms 
in this representation with the use of the Hubbard-Stratonovich transformation we find the action in the standard 
form (see, e.g., |pT[| ) 

Z[h\h]= J D[zl,z m }DA^DXexp{-S[zl,z m ,A^\}) (4) 



1 



1/T 



Slz^, z m , A/j,, A] = 2p s J dr J d d r [\(d^ - iA^)z 
o 

+i\(\z\ 2 - I) + z + hz^] 



(5) 



To describe the magnetic long-range order we have to take into account the spontaneous symmetry breaking. We 
pick up the average part that is supposed to exist for the field zn: Zi m = z m + z im , z m — (zn)S m ^. We consider 
the renormalized classical region, which is sufficient to calculate critical exponents. Then only the zeroth component 
in the Matsubara sums should be retained. Note that the below results for critical exponents are valid also in the 
ferromagnetic case since only dynamic part of the action (|5|) is different for ferro- and antiferromagnet. 
Following to [[flj we integrate out the field z from (|^) , 



Z[h\h] = J DAf,D\exp(NS ef f) 



(6) 



The effective action takes the form 



'e// 



d d k 



d d q 

gj (2nY J (2^) 
2 



4^ (z^Ro k + (k 2 S kq + R kq ) 1 {zR q0 + h q ) 



+ - 



fJ 



i\k=o ~ I z \ 2 R 0Q - (z^h k =o + zh\ Q ) - In det [kiS kl k 2 + Rk x k 2 

J ki,k 2 



with g = N/p s , 



Rk x k 2 = (fci + k2) fl Ak 1 -k 2 ,fj. + i^k!-k 2 + ^ A q ^k uf j,Ak 2 - q 

<? 

The constraint condition (0) takes the form 



(7) 



(8) 



V — 



d 2 z 



,, ia Ohq m dhq Tn 



(9) 



h=0 



where h qm = 2p s h qm . According to the structure of the transformation (|l]), the non- uniform longitudinal susceptibility 
for N = 2 is expressed as a quartic form in the fields z. For an arbitrary N we adopt the following definition: 



pkm dhl_ q m dhl +q m dh pm dh kr . 



(10) 



h=0 



Expanding (^J) in powers of 1/N yields a series for Z. Below the magnetic ordering temperature Tm the Bose 
condensation takes place and the constraint (^) in the mean-field approximation [N — > oo) [j5|6| takes the form 



gT f d d k 1 
n B = 1 - — 



{2ir) d k 2 

where the number of condensed bosons tib = \ ~z 2 plays the role of the long-range order parameter in the mean- field 
approximation. The susceptibility reads 



,„ ^ 1 2 i 1 [ d d k 1 
X (?) = ^n B S q0 + ^ J (27r)dfe2(fc + 



9 )S 



(11) 



In the next order in 1 /N we have to introduce inverse propagator for the field A 

, 4nu 



/" d d /c 1 

nA(9) ^ T i (27T)^(fc + 



5) 2 



39 



TK d A d , 4n B 



59 



(12) 



and spatial part of inverse propagator for the field A 



2 



d d k 1 /■ d d fc (2fc + g) i (2fc + g) j 



(27r) d fc 2 7 (27r) d fc 2 (fc + g ) 2 
(Mi = ^A{q)(Sij j-) (13) 



n A (g) 



4n B 



(14) 



Here 



r(d/2)r(2 - d/2)T 2 {d/2 - 1) 
2I\d- 2) 



A7 1 = 2 d - 1 7r d/2 r(d/2) 



(15) 



Note that, unlike the two-dimensional case Jill; the polarization operator II A has no "dangerous" dependence q 2 at 
small q and does not lead to infrared divergences. 

Diagrams describing 1 /TV-corrections to the one-particle Green functions and non-uniform susceptibility are repre- 
sented in Fig 1 and Fig. 2 respectively. Corresponding analytical form of the constraint equation (0) is 



1 = 



gT f d d k 1 



{2ttY k 2 

gT 



n B 1 1 - 7777 



2 

d d k 



d d k 1 
(27r) d F 
1 



[E A (fc)-S A (fc)] 
E A (fc)-£ A (fc) 



2AT 7 (27r) d n A (fc) 



fe^O 



(16) 



Here 



£ A (fc) 



^ A {k) 



T 

N 



d d q 



1 



(27r) rf n A (g) 
4-d fc 2 



1 



d(d - 2) TV A d 



(fc + g) 2 
L(k,n B ) 
1 1 



4T 



fd-1) 2 fc 2 ... , 
4^ — — ^-L(k,n B ) 



(2ir) d (k + q) 2 n A (q) 



1 

fc 2 



fc 2 - 



d(d - 2) NA d 

are the self-energy corrections (see Fig.l) calculated to logarithmic accuracy and 



(17) 



(18) 



L{k, n B ) 



ln(A d - 2 /n B ) k d - 2 < n B 
ln(A d - 2 /k d - 2 ) k d - 2 -»n B 



A being a cutoff parameter. Delta- like contribution to the susceptibility is obtained by picking out the terms with 
two condensate lines in the diagrams. They read 
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1 



d d k 1 1 



N J (2ir) d fc 4 n A (fc) 
E A (fc) -£a(*0" 







„0 



(19) 



It is natural to define the (staggered) magnetization as a square root of the coefficient at the delta-function in 
(|l9|). Note that to first order in 1/N the magnetization does not already coincide with density of condensed bosons. 
Evaluating integrals in (jl^) and ( fl9| ) yields corrections to the critical exponent (3 due to presence of singular terms of 
order of In n B . The result is 



(3=1 



,d 2 -d + 2 
! NA d 



(20) 
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Corrections to the staggered susceptibility are determined by diagrams in Fig. 2. Picking out logarithmically divergent 
in q terms at T = Tm (tib = 0) an d q — > leads to the result 

,= (d -2)(i-jy (2D 

Using the scaling hypothesis we can calculate other critical exponents. The results to first order in 1/N are summarized 
in the Table. For d — 2 + e we write down two first non- vanishing terms of the expansion in e. For comparison, critical 
exponents to first order in 1/N in the O(N) model p|Jl4| and results of the renormalization group (RG) analysis at 
d = 2 + e in the nonlinear sigma model p5| are presented too. 



Tablet. Critical exponents in different models 





d= 2 + e 


d= 3 




0{N + 1) 




RG 


0(N) 


cp n-i 


a 


-2/8+1 + 2/N 


-2/e + l + 4/N 


-2/e + 3 


-1 + 32/tt*N 


-1 + 1U/tt z N 





1 - 2e/N 


1 - 4e/N 


1 - 2e 


1/2 


1 - U/tt' a N 


V 


e-2e7(iV -1) 


e - Ae'/N 


e-2e z 


8/3n z N 


1 - 32/tt' a N 


vl/e- 1/{N - 1) 


l/e-2/N 




1 - 32/3n z N 


1 - 48/tt' z N 



The value of the exponent v for the CP model was obtained earlier from the one-particle excitation spectrum 
[|i6|,pp|] . However, only this exponent can be calculated in such a way. To find a second exponent and thereby all the 
other ones, an above-performed explicit calculation of two-particle Green functions is necessary. 

One can see that for d = 2 + e, N = 2 our results in the CP N ^ model coincide with the corresponding RG 

values. The results in the 0(3) model differ from those in the CP 1 model. This is due to that the true parameter of 

l/TV-expansion in the 0{N) model is in fact 1/{N — 2), so that higher order terms are needed to transform powers 

of 1/N into 1/(N — 2). A similar situation takes place for d = 2 where higher-order terms are needed to correct 

preexponential factor in the correlation length within the 1/N expansion for the O(N) model |fl7|| , but not for the 
CpJV-i modcl g_ 

At d = 3 the leading 1/N— corrections to critical exponents are too large for N — 2 , so that we obtain negative 
values of critical exponents. Apparently, high-order terms are required to eliminate this difficulty. Thus the mean 
field approximation which describes satisfactorily low-dimensional systems, is not too convenient starting point 
for 3d magnetic systems. At the same time, in the 0(A^)-model the 1/N corrections are small for N — 3 and yield 
physically reasonable results (see [H). 

At d — 4 — e, the propagators IT , H~l contain the small factor 1/Ad ~ e. Therefore the corrections to the critical 
exponents (^), (|2l]) are small. At the same time, the values of mean-field critical exponents (/3 = 1, r) = 2) are in 
drastic discrepancy with their values j3 = 1/2 — O(e), r\ = 0(e 2 ) calculated by renormalization group approach in the 
(f) 4 model [[l4|]. Taking into account several terms in 1/N— expansion cannot change the situation in this case. Thus 
summation of an infinite sequence of diagrams is necessary to correct critical exponents at d = 4. The situation for 
d — 4 — e is also complicated and not quite clear. It was shown by using the e— expansion Jl6|,p[ that to first order in 
e the second-order phase transition in the CP N ^ model exists in a standard sense only at 2N > 365, and at smaller 
TV the critical exponents turn out to be non-real. This contradicts to the fact of existence of a second-order transition 
in the Heisenberg model (N — 2). 

To conclude, we have calculated critical exponents in the CP W_1 modcl to first order in 1/N starting from the 
mean- field approximation |i| . For d = 2 + e the results agree with the scaling analysis |l5| and the corrections are 
small. This confirms the fact that the mean field approximation works well for low-dimensional magnets. For d = 4 — e, 
our approach, which is based on the 1/iV expansion, turns out to be poor. In the intermediate case d = 3 one has to 
investigate next-order corrections to obtain physically satisfactory results. 
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Fig.l. 1/N corrections to the one-particle Green's function of the field z. Solid line is the bare Green's function, 
dashed line is the propagator IT -1 of the field A, and wavy line is the propagator II ~^ of the field A 
Fig. 2. 1/N corrections to the longitudinal non- uniform susceptibility. 
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